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Abstract. We present a general formulation of the theory for a non-minimally coupled perfect 
fluid in which both conformal and disformal couplings are present. We discuss how such non- 
minimal coupling is compatible with the assumptions of a perfect fluid and derive both the 
Einstein and the fluid equations for such model. We found that, while the Euler equation is 
significantly modified with the introduction of an extra force related to the local gradients 
of the cnrvature, the continuity equation is unaltered, thus allowing for the definition of 
conserved quantities along the fluid flow. As an application to cosmology and astrophysics 
we compute the effects of the non-minimal coupling on a Friedmann-Lemaitre-Robertson- 
Walker metric at both background and linear perturbation level and on the Newtonian limit 
of our theory. 
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1 Introduction 

There is a general consensus that modern cosmology has entered a golden age, where in¬ 
creasingly accurate observational tools are allowing us to push our knowledge of the Universe 
beyond the limits of what was imaginable just a few decades ago. This grand project is rest¬ 
ing on few theoretical pillars: The Copernican principle, the standard theory of gravitation 
i.e. General Relativity (GR) and the big bang model of the Universe (including an inflationary 
phase). The combination of these operative frameworks has lead to a compelling coherent 
theory of the Universe which is often referred to as the concordance model or standard 
model of cosmology (ACDM) which so far seem to explain very well current observations [1]. 
Nonetheless, in the last decade there has been a growing feeling in the cosmologist community 
that alternatives to or modification of this standard model are required. 
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There are several reasons one can adduce in support of this general trend, however, there 
is a general consensus that the main issue with the current paradigm is the requirement of 
an overwhelmingly dominating dark sector in the energy/matter composition of the universe 
(about 95% of the total). Not only we do not have at the moment a definitive explanation 
for the origin of dark matter and dark energy but current observations are also questioning 
our capability of effectively describe these cosmological components with standard methods, 
e.g. by modelling them on large scales as perfect fluids. 

Indeed, approaches that attempt to describe Dark Energy (DE) as a fluid beyond the 
perfect fluid assumption [2-5] have been proposed, as well as fluid formulations of DE models 
based on scalar field [6, 7] or assuming exotic equations of state [8-10] or couplings between 
DE and perfect fluid matter [11]. Similarly, in the context of Dark Matter (DM), described 
in the ACDM model as a pressure-less dust (which can also be thought as a perfect fluid 
with negligible pressure, although the exact limit of collisions matter strictly speaking does 
not admit a fluid description), several alternatives have been proposed, mainly motivated 
by observational challenges that such standard candidate is currently facing at astrophysical 
scales [12]. In particular, alternatives to Cold DM such as Warm DM [13-15], interacting DM 
[16, 17], ultralight scalar DM or Bose-Einstein condensate (BEC) DM [18-20], have gained 
a remarkable attention thanks to their potential ability of curing small scales issues while 
remaining compatible with the ACDM on large scales (See [21] for a review of DM candidates). 
Einally, one cannot omit to say that the necessity to introduce such dark components in the 
Universe energy content seems rooted in the application of GR at scales well beyond those 
currently test with high precision [22, 23]. This has stimulated a novel activity in exploring 
extended theories of gravitation which often inclndes extra gravitational degrees of freedom 
beyond the metric (a prototype of these theories being scalar-tensor theories). 

In this paper we focus on a different approach to extend the standard paradigm. In fact, 
while all the previous alternatives rely on a modification of the dark component characteris¬ 
tics, or of GR, here we focus on a generalization of the way the former interacts with gravity. 
In particular, we will consider a perfect fluid, embedded in a curved space-time, that is non- 
minimally coupled (NMC) to curvature. In doing so we shall show that this turns out not 
only in a modification of the fluid dynamics bnt (as expected) in a subtle modification of the 
gravitational one. It is perhaps important to stress that such a modification of GR dynamics 
will be associated to a “coarse grained” regime (the introduction of cosmological perfect flu¬ 
ids and their effective coupling to gravity) rather than to some fundamental modification of 
gravitational dynamics. As snch, it does not rest on the introduction of fundamental degrees 
of freedom. In this sense this framework is quite different from proposed modifications of the 
gravitational dynamics aimed at completely avoiding the introduction of dark components. 

The use of fluids in the GR is probably as old as Einstein’s theory itself. Their simple and 
intuitive formalization has made such description of the gravitating matter a fnndamental tool 
in the understanding of the gravitational dynamics and its fonndations. Moreover, the theory 
of flnids has beneficed from its applicability to many different scales and systems [24] thus 
becoming a very well developed framework. However, even in flat space-times, a definitive 
formulation of a fully relativistic theory of fluids beyond the simplest assumptions is still 
missing. When one moves to the context of cosmology, and hence of GR, the game becomes 
even more complex as a conseqnence of the interplay with the gravitational dynamics. Indeed, 
it is still an open question how to extend the Lagrangian formulation of flat space-time 
theories when moving to curved ones. 

Non minimal coupled fluids has been investigated in the attempt of reconciling MOND 
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with DM [25], in the context of extensions to the dark matter paradigm in [26] and their appli¬ 
cation to cosmology [27]. Also the above mentioned proposals of DM condensates stimulated 
the studies of non-minimally coupled relativistic BEC [28].^ 

But why a cosmological perfect fluid should be non-minimally coupled to gravity? At 
the foundation of the fluid paradigm in GR, there is the assumption that the continuous limit 
(i.e. the limit from particle to fluid) is reached on scales small enough to be well approximated 
by flat space time. However, one may wonder what happens if the scale over which the fluid 
can be defined (i.e. the distance that a particle has to travel in order to have enough collisions 
to be thermal) is of the same size of the characteristic length over which curvature changes. 
This is very likely the case for DM, which, being almost non-interacting, has a very large free 
streaming length so that its mean free path can be of the order of the Hubble scale. More 
precisely we can estimate it from the formula 

^mfp =-^ io3 Gpc, (1.1) 

O’dm-dm Pdm 


where we have taken Udm-dm/AIdm ~ 7cm^/g [31], while Pdm is the background DM density. 
With such large mean free path it seems almost natural to include NMC terms in the action. 

A similar reasoning applies to the models based on fluids which have an intrinsic length 
scale, like the healing length of the BEC DM candidates. Indeed, these models require some 
form of self-interaction for the dark matter (keeping interactions with standard models field 
absent or negligible). This might lead to astrophysical scale condensates with potentially 
very large coherence (healing) lengths which would make dimensionally viable non-minimal 
coupling terms of the Horndeski type (see [28] for a more detailed discussion). 

Besides these heuristic arguments and their application to the DM problems, the idea 
of a non-minimally coupled fluid is interesting per si as it provides further insights on the 
dynamics of fluids in GR. In fact, this kind of models show several intriguing and distinctive 
feature that make the study of their phenomenology interesting. As we shall see a NMC fluid 
sources the gravitational potential also via the gradients of its density distribution. Moreover, 
the equations of motion for the fluid show interesting features. In fact, the NMC introduces 
an extra force in the Euler equation that is related to the gradients of the curvature, while 
the continuity equation is modified in such a way that it is still possible to define conserved 
quantities along the flow lines of the fluid. All these modifications have clear relations with 
observables and, especially in cosmology, can be easily compared with data. 

The paper is organized as follows: in section 2 we review the theory for a perfect fluid in 
curved space times while in section 3 we extend these results to the case in which the fluid is 
non-minimally coupled to curvature. In particular we derive the Einstein equations and the 
fluid equations for both a conformal and a disformal coupling and discuss their properties. 
In section 4 we analyse the cosmological consequences of the model at both background and 
linear perturbation level while in section 5 we derive the Newtonian limit of the NMC fluid. 
In section 6 we discuss the possibility to have more than one fluid which is NMC and the 
effects that such coupling can have on the DM-DE interaction. Einally, in section 7 we draw 
our conclusions. 


^The idea of NMC matter has been also used as an alternative to DM itself [29, 30]. 
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2 The action for the relativistic perfect fluid 


A fluid is defined to be perfect when it has isotropic pressure, no viscosity, and (in relativity) 
no heat conduction in the comoving reference frame. As such it can be described by a set of 
five independent variables, namely the three components of the fluid four velocity and by 
any two scalar functions describing the thermodynamic state of the fluid [24, 32]. However, it 
is not possible to construct a valid variational principle only with the previous functions unless 
this variation is constrained [33]. One possible way to constrain the variation is to introduce 
Lagrangian multipliers in such a way that the correct fluid equations can be recovered once the 
properties imposed by those constraints are applied [34]. Our action will be then constructed 
from a scalar F, function of particle number density n and entropy density per particle s, i.e. 
by an equation of state, and from a set of Lagrangian multipliers. The choice of the latter is 
based on a minimal constrain approach: The variation is constrained in such a way to give 
the correct equations of motion for a perfect fluid. Practically speaking, this amounts to 
require the particle number conservation and the absence of entropy exchanges between flow 
lines. It will turn out that a more effective description of the action can be done in terms of 
a vector density field dehned as 

= ( 2 . 1 ) 


which in particular gives n = \ J\/^/—g and = J^/\J\, being respectively a scalar and a 
pure vector and with j Jj = ■ The action for the perfect fluid is then [34] 


^fluid = J d^x^F{n, s) + + sV^e + . (2.2) 


The first term is the equation of state of the fluid while the second one it contains the required 
constraints. The Lagrangian multipliers, (p, 9 plays the role of thermodynamic potentials 
and are associated with the chemical free energy and the temperature respectively. The last 
constraint is instead required to ensure that the flow lines do not cross, i.e. there is a one-to- 
one mapping between the Lagrangian and the Eulerian coordinates at any two fixed spatial 
hyper-surfaces. In particular, the and A = 1,2,3, are the Lagrangian coordinates of 
the fluid. This two part structure has the particular advantage that it divides the action into 
a metric dependent and metric independent part. This help in distinguishing the properties 
of the embedding of the fluid in a curved background from the symmetries that must be 
enforced independently of the geometry. 

In fact, the variation with respect to the metric gives a stress energy tensor (SET) for 
the fluid 


rril-LU _ 

^ 7 ^ dg, 


SSnuid ^ p p 

on 


fll/ 


(2.3) 


This has the form of a perfect fluid SET and defining p = T^'^u^Uy and p = it can 

be cast in the standard form = {p + p)u^u'^ + pg^^'^ with the identifications 


^ dF ^ 

p= -F , p = n— - F . 


(2.4) 


This also tell us that the Lagrangian for the perfect fluid is the energy density measured by 
an observer comoving with the fluid. ^ Adding the standard GR kinetic term for the metric 

^Equivalent Lagrangian densities can be obtained adding surface terms and using the thermodynamic 
relations for the fluid. 
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and varying it with respect to the metric gives the Einstein equation 


= (2.5) 

The equation of motion for the perfect fluid can be obtained by taking the covariant derivative 
of the Einstein equations. With this procedure we find the standard continuity and Euler 
equations 


p + e{p + p) = Q, ( 2 . 6 ) 

{p + p)u^ + hW'^p = Q. (2.7) 

where we have defined () = tt^V^() and 9 = It is interesting to note that this is a set of 

four equations for five variables. Usually, this issue is solved by assuming an equation of state 
for the fluid which relates the pressure to the density. However, recalling that p = p{n,s), 
the continuity equation can be expanded as 

+ + = p.8) 

We now see that there is one missing equation to close the system. This is provided by the 
request that the particle number density is conserved along the flow. This amounts to say 
that the first parenthesis is independently zero so that the conservation of entropy follows 
automatically and the system is closed [35]. 

Until now, no use of the constraints has been done. However, we will discuss now how 
their role is fundamental in order to obtain the correct equations of motion for the fluid 
when those are derived by varying the action with respect to the fluid variables. In fact, such 


variation leads to the following set of equations 

<5 S 

0 = = pUfj,+ \/+ s\7^9 +/3 a^ fiO:^ , (2.9) 

0 = ^ = -V^J'^, (2.10) 

0 = ^ =-V^(sJ^), (2.11) 

0 = ^ = -V^(/3aJ^) , (2.13) 


According to the dehnition of the continuity equation is obtained by combining together 
the equation for the entropy (2.11) and number density conservation (2.10) while the Eu¬ 
ler equation is obtained by differentiating equation (2.9) with respect to Vj/, then anti¬ 
symmetrizing it and projecting it along the mixed direction [34]. 

The fluid equations can be rephrased in a more thermodynamic form contracting the 
first equation with the four velocity and using the relation between and the other fluid 
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variables (2.1) 


fi = f + Ts, (2.15) 

V^(nu'^)=0, (2.16) 

= (2.17) 

u>^V^Pa = 0, (2.19) 

= 0 , ( 2 . 20 ) 


where = (p + p)/n, f = {p + p)/n — Ts. In this form the equations give the relation 
between enthalpy p, temperature T and free chemical energy / and express the conservation 
of number density and of entropy while the last two equation are more of geometrical nature 
and tells us that the Lagrangian coordinates are constant along the flow lines. 

As a final remark we discuss in more detail the role of the two parts of the fluid La¬ 
grangian. The fluid being perfect is related to the fact that in its metric dependent part no 
derivatives of the fluid variables appear. Hence, its SET can only be of the perfect fluid form. 
The role of the constraints is to enforce the properties of the fluid we know from the conser¬ 
vation of the SET when doing the variation with respect to the fluid variables. Of course, 
the two parts are related in order to give a consistent theory so, for example, inclusion of 
diffusion terms will need a change in both parts of the Lagrangian. Hence, it is not trivial 
how to generalize this action to more complex situations. We will see in the next paragraphs 
how one can obtain a more general action for a fluid still compatible with the constraints. 

3 Non-minimally coupled fluid 

As we have discussed in the introduction, the minimally coupled perfect fluid paradigm may 
not be well motivated in GR for fluids with mean free paths of the same order of the curvature 
scale and we have pointed out how an effective coupling between fluid variables and curvature 
may rise. 

However, when moving from flat to curved space-times, there is no general prescription 
on the kind of couplings between curvature and fluid variables and their derivatives so that an 
a priori selection criterion must be adopted. Here we take advantage of the natural derivative 
structure of fluids to argue that there are only two non-minimal couplings that can be added 
to the perfect fluid Lagrangian without introducing derivatives of the fluid variables;^ a 
conformal coupling, Fc{n,s)R, and a disformal one, F^^n, 

Hence, the most general action that contains non-minimal couplings and is still zero 
order in the fluid derivatives is [26] 

I- 

S = / d'^Xy/^lR + acFc{n, s)R + a^Fain, s)R^^u>"u‘'] + 5fiuid , (3.1) 

^Modulo surface terms. Notice also that by derivatives of the fluid we mean terms that cannot be eliminated 
via integration by parts. 

■^One may argue that, as curvature terms contain derivatives of the metric, such NMC terms may be 
thought as fluid derivatives contributions hence requiring an extension of the model to include all possible 
derivative terms up to second order. However, those terms manifestly break the perfect fluid assumption and 
will for this reason not included. 
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where S'fjuid is given by (2.2). Notice that the associated total Lagrangian can be rewritten 
as 


C = + F(n, s) 


1 + 


Mil 


Fc{n,s) Fd{n,s) 

Oic.^r, - rR + - -Rfj^yU^U 


F{n, s) 


F{n, 


(3.2) 


which can be seen as a shift in the equation of state of the minimally coupled fluid. In the 
simplest case oi Fc = F and aa = 0 the shift is constant from the fluid variable point of view 
and hence we expect to see changes in the equation of motion for the fluid to appear as shifts 
and rescaling (See appendix A.l) . Of course, the Einstein equations will be modified more 
signihcantly. In this sense, there will be differences in the conserved fluid quantities and the 
ones sourcing the gravitational potentials. However, as we shall discuss in more details below 
which one to choose will be mostly a matter of convenience. 

Before moving the the details we note that the presence of a non-minimal coupling in 
the action will produce a coupling between second derivatives of the metric and the fluid in 
the Einstein equations so that there is no unique way to dehne, from the metric equations, 
the SET components as seen from an observer. In fact, the Einstein equations will take the 
general form 

M|i(V^)G^, = (i?, rA, d^lj, dd ^), (3.3) 

where '0 stands for fluid variables and where contains both the minimally coupled perfect 
fluid SET and the corrections due to the NMC. One may then envisage at least three different 
schemes for writing these equations: 


• Define an effective Planck mass M| and an effective SET; 


• Eix the Planck mass to be constant and dehne an effective SET that will contain the 
Einstein tensor; 


• Use contractions of the Einstein equations with the metric and the fluid four velocity to 
get rid of the curvature on the right hand side. However, the presence of second order 
derivatives of fluid variables will still mix gravitational and huid degrees of freedom. 


Of course there is no physics into this and the choice of which prescription to use will 
be mostly determined by the observables chosen to compare the model with data or by the 
numerical procedure used to integrate the equations.^ Here, we will pick the hrst possibility 
which was recently used in [36] so that the Einstein equations are 

= (3.4) 

while the fluid equations are derived by taking the covariant derivative of such equation, 

. (3.5) 


Notice that, being in the Jordan frame, any other matter species which is minimally coupled 
will be separately conserved [37]. 

Einally, as a consistency check, in appendix A.l we have also computed the fluid equa¬ 
tions by varying the action with respect to the fluid variables, including the perfect fluid 
constraints, and found the same equations obtained from the conservation of the Einstein 
equations. 

physical definition of the effective Planck mass can be given by looking at the coefficient of the graviton 
kinetic term appearing in front of the second derivative term in the expansion of the action for perturbations 
around some background. 
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3.1 Conformally coupled perfect fluids 

As a first investigation we consider the case of a perfect fluid which is coupled to the Ricci 
scalar. In this case the action is 

]\,t 2 n 

Sc =—^ J d^xV^[l + acAc(n,s)]R + Sfluid, (3.6) 

where Aguid is given in (2.2) and Fc{n, s) is a new arbitrary function of the fundamental fluid 
variables. 

The variation with respect to the metric gives the Einstein equations for the conformally 
coupled fluid 

( dF OF \ 

+ etc 2^ j {gni^OFc — V^V uFc) = 0 . 

(3.7) 

The contributions coming from the NMC appear as a modification to the Planck mass that 
is now fluid dependent, as a correction to the perfect fluid structure of the SET and as a 
derivative contribution. At this point one may want to follow the same strategy used for 
the minimally coupled case in order to identify the thermodynamic quantities of the fluid, 
like energy density and pressure. However, now there is no general prescription that allows 
one to unequivocally dehne the fluid variables from the Einstein equations. According to the 
discussion presented above, we then define an effective Planck mass 

= APp[(l + OicFc ), (3-8) 


and an effective SET 


^eff 


QfiiyF h^i/Tl I “1“ 


AP|i ^dFr 


R 


dn 


o^cAfpj {gfj_i/\Z\Fc — fj.'^i/Fc) , 


(3.9) 


so that the Einstein equations are 


= r; 


-eff 

^ IJ.U ■. 


(3.10) 


where and T®® are given in (3.8) and (3.9) respectively. 

The fluid equations can be derived by taking the covariant derivative of the Einstein 
held equations. After some algebra we get 


VjyE+ac 


Ml 


PI 


RVuFr-n 


dF M^.dFA,. ^ ^ { dF dFr 




= 0 . 


(3.11) 


Projecting this equation along vy gives 


dF 




.dF, 


]in + 9n)+( — 


dn 


dF 


.dFr 


ds 




s = 0, 


(3.12) 


As for the case of the minimally coupled perfect fluid, we need to impose an extra equation in 
order to close the system. The structure of equation (3.12) and its similarity to the analogous 
one for the case of a minimally coupled fluid suggest that we take again as extra equation the 
conservation of the particle number density along the flow. Such imposition is also the one 
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that makes compatible the above system of equations with the one derived in appendix A.l.® 
As a consequence of this, any scalar function of n and s will obey the following equation 

dT 

T + en— = 0, (3.13) 

on 

with = J^(n, s). In analogy with the minimally coupled perfect fluid this allows us to write 
a formal continuity equation defining a density pjr = and a pressure pjr = —ndJ-jdn + J- 
so that 

pjr + 9{pjr + pjr) = 0 . (3.14) 

This is interesting as it implies that even with the NMC it is possible to identify quantities 
that are conserved along the flow line.^ In particular, F{n,s) and Fc{n,s) both (or any 
combination of them) are conserved and hence can be used to define fluid variables. 

The projection onto the orthogonal hypersurface gives 


n 


dF 

dn 


+ cx, 


2 dn 


u — 




^ dF 

F - n- - ac 

On 


M|i 


OF, 

dn 


n-FAR]+ac^Fch‘^’'V^R = 0. (3.15) 


We can see that the NMC has the effect of modifying the weight function suppressing the 
spatial gradients sourcing the time derivative of the fluid velocity, the pressure contribution 
and add a new potential related to the local value of the Ricci scalar. Schematically we can 
rewrite the previous equations in a form that is closer to the one for a minimally coupled 
fluid (2.7), namely 

{ptot + ptot) (ptot) - ac^p.h^'^V.R = 0 , (3.16) 

where Pc = —Fc while ptot and ptot are collective labels for the various coefficients appearing 
in equation (3.15) and read 


Ptot — F + —MpiRFc 


Ptot = n- 


. dptot 
dn 


Ptot • 


(3.17) 

(3.18) 


This is interesting as it shows that even in the simple case of conformally coupled fluids, 
the NMC has a clear and distinguishable feature: It introduces an extra force term. In fact, 
while the pressure renormalization can be in principle mimicked with other exotic fluid, the 
presence of a new force term related to the local value of the Ricci scalar, is a signature of this 
model. By means of the trace of the Einstein equations and taking terms up to 0(ac), we 
can relate this term to the gradients of trace of the SET of all the matter helds. Hence, this 
force will be effective not only in the presence of gradients of the NMC fluid but whenever 
any matter distribution is sufficiently inhomogeneous. 

®This is not a completely arbitrary choice. Indeed, the coupling between the variables n and s and curvature 
is ‘scalar’ in the sense that these variables are coupled to scalar quantities, namely R and (R) and hence we 
do not expect in the equation of motion terms able to generate a spatial particle flow. The situation would 
be different if we had considered a NMC of the form Rij,^\7i^F{n, s)V'F{n, s) which has a clear direction 
dependence. 

^Interestingly, an analogous result has been found in the case of perfect fluid coupled to a scalar field [11]. 
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3.2 Disformally coupled perfect fluids 


We now consider the case in which the NMC is of the form Fd{n, s)R^i,u^u'^. The action for 
this model is 


Sd 



2 


{R + adFd{n,s){R)) 


T ‘S'fluid ) 


(3.19) 


where iSfluid is given in (2.2) and Fd{n,s) is a new arbitrary function of the fundamental 
fluid variables, while (R) = R^yU^^F^. Even in this case no derivatives of the fluid variables 
are present in the action. Here, however, one could re-express the NMC coupling using the 
commutation rule of the covariant derivative 


R^.yU^^F' = (V^V^u^ - , (3.20) 

so that an apparent higher order derivative term has appeared. However, given that this 
combination is zero in flat space time it is more reasonable to consider such interaction as a 
gravity-fluid coupling rather than a thermodynamic property of the fluid. 

Taking the metric variation one gets the following Einstein equations for the disformal 

fluid 


dF M‘^ / 

FtpiG^1/ F “h 'il^lly^Tl^— ^ ^ Fdi^lVji^g^y + ‘lU^Uy) 

2V„V(^^“ -h + g^u{R)n^F + = o, (3.21) 


where the round brackets stand for symmetrized of the indices and where we have introduces 
for simplicity the tensor 

F’' = Fduy ’'. (3.22) 

Notice that in this case the Planck mass is constant so that in principle we could dehne 
the observed fluid variables by decomposing the Einstein equations in its effective variables 
but again this would hide in the matter variables curvature terms that will ultimately mix 
gravitational and fluid degrees of freedom. Eor this reason and to be consistent with the 
formulation of the conformally coupled case we stick with our notation. We then have a 
constant Planck mass = Mpi and an effective SET 


dF^ i\'(i ^ f 

^ puV F 'U^lly^Tl-^— Old ^ ^ F(j^(^Fi}j {^g^y -\- ‘lU^Uy^ 

+ - Ut^y + g^yVa'^gt^'^ - gi,y{R)nFF^ _ {R)ui,UyF^'^ 


Ml 


dn 


dn 


,(3.23) 


so that the Einstein equation for the disformally couple fluid are 

M^,G^y=Tf,. 


(3.24) 


where the effective SET is given by the expression above. We now move to the computation 
of the fluid equations of motion. Again we take the covariant derivative of the Einstein 
equations and project the result along the direction parallel to the fluid four velocity and 
onto its orthogonal hypersurface. Quite remarkably for the hrst one we get 


Mgi dFd 
\dn ^ 2 dn 



(h -|- 6n) + 


(dF , Ml, dFd 



s = 0, 


(3.25) 
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which again tells us that the particle number density and entropy per particle are conserved 
along the flow lines, with the conclusion that any scalar function of these variable obeys the 
continuity equation (3.14). Even if surprising, this result is in line with the expectation that 
our NMC action is compatible with the perfect fluid assumptions. The spatial projection of 
the fluid equations gives the modified Euler equation 


n 


dF 

dn 


- adiR) 2Fd - n 


dn 


+ = Od 


20 ( Fd-n 


if + 

dFd 


n 


OF 

dn 


-F + a0n^-F,]{R) 


dn 




(3.26) 


This equation is more complicated than that for a conformally coupled fluids, but still we can 
see some common features. There is a modified coefficient in front of the time derivative of 
the fluid four velocity, a modified pressure term and an extra force related to the Ricci tensor. 
The new feature of this case is that on the right hand side of the equation we have two terms 
that act as source for the equation. These terms represent a new kind of contribution which 
is not related to spatial gradients and hence represents another novel feature of the NMC. 
Notice also that any component of the acceleration will in general receive contributions 
from all the others via the contraction h^^Rp^if present in the last term of the second line. 

Finally, we note that the last term in the second line, despite the appearance of a 
derivative of the Ricci tensor along the flow line of matter, is not introducing any new degree 
of freedom. In fact, it can be rewritten as 

(r^ 0U^^ = - V0V^uf^^ , (3.27) 

which shows that of the three derivatives at most two are time derivatives so that only second 
time derivatives of the metric will appear. 


3.3 Count of degrees of freedom 

In the previous sections we have seen how the coupling of fluid variables to curvature terms 
implies the appearance of second derivatives of such variables in the Einstein equations. 
Hence, an important point to address in this class of models is whether the NMC introduces 
unwanted degrees of freedom. Stated in another way, one needs to check that the standard 
relativistic hydrodynamics is not spoiled when introducing such non-minimal couplings. 

It is the aim of this section to discuss this point for the two cases of conformally and 
disformally coupled perfect fluid, providing arguments that point towards the absence of such 
extra degrees of freedom. 


3.3.1 Conformally coupled perfect fluid 

In the case of conformally coupled perfect fluids the answer is straightforward. The equations 
for this case are (3.7), (3.15), (A.2) and (A.3). In the Einstein equations, second derivatives 
of the function Fc appear which may lead to unwanted second time derivatives of the fluid 
variables n and s. However, the conservation of the particle number density and of entropy, 
the last two equations of the system, implies that we can always eliminate such derivatives 
and hence no extra degrees of freedom are introduced in this model. 
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3.3.2 Disformally coupled perfect fluid 

In the case of disformally coupled perfect fluids the situation is more complex. In fact, in the 
Einstein equations (3.21) appear derivative of the tensor t^i, = Fd{n, s)Uf^Uu- When those 
are expanded they will give second derivatives of the scalar functions n and s as well as of 
the fluid four velocity. The former can be eliminated by the same arguments used for the 
case of a conformally coupled fluids as also in this case conservation equations holds for those 
variables as can be seen from equations (A.13) and (A.14). However, such constraint does 
not exist for the fluid four velocity. It is hence important to check that the Einstein equations 
are actually free of second time derivatives of the fluid velocity. In order to show this we 
will make use of the natural space-time foliation provided by the vector field and project the 
Einstein equations accordingly.® 


Time projection Projecting the Einstein equations along the tensor it gives the 

following scalar equation 


A4^{G)+F-^M^^Fd{R) + ^M^^ =0 (3.28) 


Od 


We now proceed with the analysis of the last parenthesis which is the one containing the 
second derivatives of the fluid variables. As a consequence of the unit norm of the four 
velocity, the following relations holds 


U^VaUfj, = 0 , = -V , (3.29) 

and hence, by expanding the derivatives in the last parenthesis and keeping only second 
derivatives of the velocity, we get 

- 2u^u^V„V^t^" + ~ 2 {iFf + 2 {Vo^upf + (R) . (3.30) 


Hence, no second derivatives of the four velocity are present. 


Space-time projection We now project the Einstein equations along the mixed direction 
where h is the spatial projector onto the hypersurface orthogonal to the fluid velocity. 
This gives 

^ = 0 . (3.31) 

Again, considering only the derivative term and keeping only the second derivatives we get 






(3.32) 


which, again, does not contain any second time derivative of the four velocity. 


Spatial projection This time we project the Einstein equation onto the orthogonal sub¬ 
space. This amount to contract the Einstein equations with the tensor . The result 

is 


dF\ 


+ [-F + n^)h^P + ^MMR) 


dn J 


-Fd + 


dn 




+ [-2V„V(^t,)“ + □ V + = 0 . (3.33) 


“This is the standard procedure to define the observed energy density, momentum and stresses of a system 
as seen by an observer whose world-line is defined by a vector u''. 
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Again, extracting only the second derivative from the last parenthesis of the above equation 
leads to the following terms 






~ -h^^hP'^u^Va + V^u^) + 2h^P6 . 

(3.34) 


The last term does indeed contain second derivatives of the fluid four velocity. In fact 
9 = However, thanks to the unit norm constraints (3.29), it is easy to show that 


9 = , 


(3.35) 


which tells us that the expansion 9 reduces to a spatial derivative of the fluid four velocity. 
Hence, the term 9 actually contains only one time derivative. We thus conclude that, as for 
the previous components, also in this case no higher time derivatives of the fluid four velocity 
are generated by the NMC. 

In this section we have made a hrst analysis of the number of degrees of freedom for a 
NMC fluid and concluded that despite the appearances no extra degrees of freedom seems 
to be generated by the coupling. Of course, only an Hamiltonian analysis can say the final 
word on this aspect but, as we will see below, our intuitions are conhrmed by the analysis of 
the linearly perturbed equations around a Friedman-Lamaitre-Robertson-Walker (FLRW) 
background. It is indeed very likely that the results derived in this section will hold also 
when a more complete analysis will be done. 


4 Cosmology of a non-minimally coupled perfect fluid 

A crucial test for all the modified gravity models comes from their ability to match the 
background evolution of the Universe and the process of structure formation. Even if a 
detailed analysis of the cosmology is beyond the scopes of this work, we here derive the main 
equations and discuss the modifications produced by the NMC. To do so we consider a FLRW 
metric dehned by the usual squared line element 

ds^ = —(1 + 2^)dt^ + a^(t)(l — 2T)dx • dx, (4.1) 

and apply it to the Einstein equations (3.7) and (3.21) respectively as well as to the fluid 
equations (3.14), (3.15) and (3.25). 


4.1 Cosmological background 


We start our investigation by looking at the equations for the background, setting the grav¬ 
itational potentials $ and T to zero. As a first step we recall that from (3.14) it can be 
inferred that the background evolution of the fluid equation of state F as well as that of 
the NMC functions Fc and F^ is determined once specific equation(s) of state are taken. In 
particular, we have 

Pi =-2>H{pi+pi), (4.2) 


with Pi = —Fi and pi = Fi — dFijdn. Recall that the NMC functions are dimensionless and 
hence their density and pressure are to be considered as divided by some reference scale. 
Einally, as it will be used later, we derive also the expression for the second time derivative 
of the density 


h = 


-3H +977^(1+ 0?) 


{PiFpi) , 


(4.3) 
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where cf = dpi/dpi is the speed of sound for the i-th species. For simplicity we further 
assume that the minimally coupled fluid has negligible pressure. In this limit, equation (2.4) 
implies F (x n so that one can expect that Fc{n, s) = Fc{F, s) as well as Ffi{n, s) = Fd{F, s). 
With the identihcation F = —p, being the latter the minimally coupled energy density, one 
can deduce that the NMC function will be generic functions of p and s. 

With these details at hand we now proceed with the investigation of the cosmological 
background evolution of our NMC model. 


4.1.1 Conformally coupled fluid 

For the conformal NMC the Einstein equations reduces to the following system 

^2 ^ 8^_ p _ 

3 1 + ac(2pc + 3pc) 


(4.4) 


H = -AttGp 


1 + acjpc - 3(pc +w)Cc) 

(^l + ac{2pc + 3pcfj 


(4.5) 


From these equations we see that the effects of the NMC reduce to a rescaling of the minimally 
coupled density p.® It is important to notice that the denominator 1 + ac{2pc + 3pc) has to 
be strictly positive in order to avoid singularities and the flip the gravitational constant sign. 
Of particular interest is the sign of the expression between brackets in the second equation. 
If the NMC terms start to dominate at late times they can flip the sign of a, thus producing 
a positive acceleration even if the minimally coupled matter species satisfy the strong energy 
condition. 

However, in order to reproduce the observed matter and radiation dominated eras, the 
NMC contributions have to become negligible as we move backwards in time. This can be 
easily achieved if, for example, the equation of state of the NMC fluid is such that the density 
Pc is a growing function of the scale factor. 


4.1.2 Disformally coupled fluid 

For the disformal NMC the Einstein equations are 


^2 ^ 8^_p_ 

3 1 - ^Udipd + 3pd) ’ 


(4.6) 


^ l+«d(Pd-i(Prf+Pd)4) 

^ (1 - \ad{pd + 3pd)Y 


(4.7) 


Similarly to what we found for the conformal case, the NMC rescales the energy density 
contribution and, in particular, we have to require that the denominator in the two equations 
is never zero. In the second equation, the denominator is always positive while the numerator 
can freely change its sign so that also in this case the NMC may change the sign of the 
acceleration. 


®If we want to include other minimally coupled matter species it is sufficient to replace p with the sum of 
all the other components. This also tells us the universal nature of the NMC as its modifications affects all 
the matter species. 
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4.2 Linear perturbations 

We now proceed with the investigation of the dynamics of linear perturbations. We define the 
density perturbation as p{t, x) = p{t) + 5p{t^ x). Thanks to the conservation of both number 
density and entropy, the perturbations in the densities p, pc and pd obey the following 
equation 

6pi + 3H + cj) 6pi - 3{pi + pi)4> + {pi + pi)V • n = 0 , (4.8) 

where c, = dpi/dpi, which is exactly the GR result as expected. 

In the next sections we will derive the remaining equations for the two cases of confor¬ 
mally and disformally coupled fluids. 

4.2.1 Conformally coupled fluid 

We now derive the exact expressions for the linear perturbation equations for a conformally 
coupled perfect fluid. We will first derive the equations for the gravitational potentials 
and T and then the one for the spatial components of the fluid velocity. 

Euler equation 


Esciu^ + Egc2H'iE -|- 3oicEsc^^Spc — 

= MPc + Pc) + 4V*^) + 3V*T + + 2V*V2$ 


(4.9) 


where the background coefficients Egd are reported in appendix B. Notice that the left hand 
side of the Euler equation has the same structure of the one derived in GR, only with modified 
background coefficients. On the contrary the right hand side is represent a characteristic of 
the NMC. It shows how the velocity is now sourced by a complex combination of derivatives 
of the metric potentials. 


Time time component of the Einstein equations 


— [1 -t- Oc (2pc + 3pc)] -1-2(1 — acPc) — K^6p = 


— Qir 


—3H‘^ [2 + 30^) 6pc + 3H{pc + Pc)'^ • u — V^dpc ■ (4:.10) 


The left hand side of this equation is formally identical to the GR expression, only with 
modified cefficients. The right hand side represent instead a novel feature of the NMG: 
gradient of the fluid variables are acting as a source for the gravitational potential. 

Space space component of the Einstein equations 


ac'^iVjSpc - (1 - acPc)(ViVj4> - VjVjT)) 


-\- h- 




(2E,i4> + 2Es2^ + 2E,3T + + 2^,5^ 


-|-(1 -|- q;cPc)V^<I? — (1 -|- otc{Pc + 2pc))V^4'] — 

= hij [+acV‘^5pc + acEsrSpc + 9acH‘^{pc +Pc)5cl 

-2acEsQV ■ u -F adpc + Pc)^ ■ u 


(4.11) 


^°This is something expected in any NMC model and has been dubbed kinetic mixing. See [38] for an 
analysis in the context of Horndeski theories. 
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where the exact expressions for the background coefficients Esi are reported in appendix B. 
From the first line is clear that anisotropic stresses are present and are related to gradients of 
the non-minimal coupling function Fc = —pc- This represents a first difference with respect 
to GR. The part proportional to the spatial projector hf^i, contains on the left hand side, 
the standard GR terms with modihed ceofficients, while the right hand side represent the 
genuine contributions coming from the NMC. 

Space time component of the Einstein equations 

p + 3ac(2i?^ + H)[pc + Pc)^ Ui + + V^T) ^1 H—^(/^c + = 

= Oc - 3Cc)Vi(5pc + (Pc+ Pc)ViV •-u)^ . (4.12) 

Also in this case we have on the left hand side the standard GR terms with modihed back¬ 
ground coefficients while on the right hand side new effect related to the gradients of the 
huid variables are present. 


4.2.2 Disformally coupled fluid 

We now move to the analysis of the linear perturbation equations for a disformally cou¬ 
pled perfect fluid. We will hrst derive the Euler equation and then the equations for the 
gravitational potentials and T . 

Euler equation 

Ediii^ + Ed2Hu^ + - 3ad(^" + H)^^5pd = 

= Od (H{pd - 3pd)V*4> - QH{pd + 3pd)V*T 

+ {pd - 3pd)V'T - {pd + (4.13) 

where the background coefficients Edi are reported in appendix B. 

The structure of this equation is similar to the one obtained in the case of conformally 
coupled fluids. The right hand side contains terms that are also present in the GR case, only 
with modified background coefficients. On the right hand side instead we And the derivative 
interaction typical of NMC theories. 


Time time component of the Einstein equations 

- 4V"T + 2t^{5p + 2p^) + 12FT(1 - ^(pd + 3pd)) = 


= ad 


H{pd - 3pd)V ■ u - V^6pd - 3^2 1 + 3c^ 6pd (4.14) 


Space space component of the Einstein equations 


ad ( -H{pd + ^Pd)'^^ + Pd 


V iUj 


ViVj4> - ViVjT 


+ hij ( Esdi^ + Esd2^ + TEsds + '^EsdA + Esdb^ + Egde 






— a^hij I Egd7 


—— + ^(Pd - Pd)'^—^ + Esdshpd + ^{Pd+Pd)Scl ) (4.15) 

a 2 a 2 


where the coeffcients Egdi are reported in appendix B. 
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Space time component of the Einstein equations 


{2k^{p + p) + adH{pd - 2>Pd) - ‘iadH‘^{pd+Pd))ui 

+ (4 — 3ad(pd + Pd)) = 


— “d 


{2Vi6pd + 3Vi6pd) - {pd + PdWN ■ u 


(4.16) 


We conclude this section with a comment on the number of degrees of freedom. As it 
has been discussed in section 3.3, the model under consideration should not propagate more 
degrees of freedom than those of standard fluid dynamics. Indeed, the investigation of the 
linear perturbations around a FLRW metric confirms that result: all second time derivatives 
of the velocity that appears from the various terms in the disformally coupled case exactly 
cancel each other. Notice that this choice for the metric is much more sensitive to extra 
degrees of freedom than maximally symmetric space-times (Minkowski or de Sitter). 


5 Newtonian Limit 


To properly discuss the Newtonian limit, it is important to carefully work out the weak field 
limit of equations (3.7) and (3.21). This was computed in [26], so here we will report only 
the main results, referring to interest reader to the cited paper. 

In this limit, the fluid equations of motion for both the conformally and the disformally 
coupled fluid formally reduce to the standard continuity and Euler equations as it can be 
easily seen by noting that any combination 72 x T) ~ (dvrG)^ with i = c, d is a post-Newtonian 
contribution to the equations. However, the dynamics of the fluid is not equal to that of a 
minimally coupled fluid as it will receive contributions coming from the modified gravitational 
potentials. 

The NMC effects remain in the gravitational equation. In particular, the Poisson equa¬ 
tion is modified, and the GR equality between the two gravitational potential is broken 
so that anisotropic stresses are expected. As we have said, any combination of the kind 
RxFi ~ (dvrG)^ is a post-Newtonian contribution to the equations and hence will be dropped. 
Then we define 

Qiiu = Viiu + Ifiu , , 7 = , (5T) 

so that the modified Einstein equations in the weak held limit in the transverse gauge for 
the conformally and disformally coupled fluid read respectively 


Mpi 

2 




-OcMpi {r]f,t,nFc - d^duFc) + g^u[F -n 


dF 

dn 


dF 

- n—u^Uy , 


-A/pi 

2 




— Q;d 






+ 9^lv 


dF\ 


dF 


n—u^u 

dn 


V ) 


(5.2) 

(5.3) 


where 

= 6l5lUFd - Sldod^Fd - SldodyFd + p^^M^Fd . (5.4) 

As one immediately sees, the effect of the non-minimal coupling is still present, even in 
the weak held limit, and the fluid is not behaving as a perfect fluid in Minkowski space-time: 
the non-minimal coupling has generated a SET which contains additional terms, constructed 
out of the derivatives of the fluid variables. 
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Putting everything together, and considering the static, non relativistic limit (i.e., the 
—>■ oo limit), we get the Poisson equation for the Newtonian gravitational field in the 
conformal and disformal case: 

(p - OcV^Fe) , (5.5) 

(p - , (5.6) 

where P) = P)/(47rG). The Newtonian potential has as sources not only the mass density p, 
but also derivative terms. In this sense it will depends not only on how much matter there 
is at one point but also on how it is distributed. 

The spatial part of the Einstein equations give 


= A-KGrjij p - OcV^p; 
= dvrGpi,- \p - 


+ SttGoc 


rjijV'^Fc + didjFc 


Making use of the Poisson equation the two spatial equations can be rewritten as 


(5.7) 

(5.8) 


+ SvrGac + didjFc 

VH,j = pij , 


(5.9) 

(5.10) 


which interestingly shows how, for a disformal coupling, there is only one gravitational po¬ 
tential, while the conformal coupling can be responsible for anisotropic stresses. 

To summarize, we have seen how both forms of NMC produce a derivative correction to 
the Poisson equation but while conformally conple fluids will excite extra gravitational degrees 
of freedom with respect to GR and will have anisotropic stresses, a disformally coupled fluid 
is described by a single gravitational potential and has no anisotropic stresses. 

A comment is also in order about the presence of the Laplacian in the Poisson equa¬ 
tions. Indeed, this implies an algebraic relation between the gravitational potential and the 
NMC functions which are ultimately related to the matter distribution. As it was noted 
in [39], this implies that any sharp change in the matter distribution can be imprinted in an 
analogous change in the gravitational potential thus being potentially dangerous for the kind 
of models under investigations. However, in our case, the NMC matter has a much shallower 
distribution as compared to standard baryonic matter investigated in the cited paper, being 
related to the description of the dark component of the Universe. Also, sufficiently smooth 
functions Fc or Fd can avoid snch issues. Hence, we do not expect that this feature of the 
modified Poisson equations will practically affect the physical viability of this class of models. 

Finally, it is perhaps worth noticing that that, when considering both conformal and 
disformal couplings, precise cancellations can happen leading to a dynamics which would 
be very similar to the GR one. For example one might consider the case in which the two 
contributions to the Poisson equation from these couplings cancel each other thus leaving 
the modihcations of gravity only to the spatial part of the potentials. Indeed this happens 
for the particular combination 


Olr — 


Q!d 


^ and Fc = Fd, 
which corresponds to the following NMC 

T^nmc = olrgFrg {—R + G^yU^vF) 


(5.11) 

(5.12) 


which, remarkably, involves exactly a coupling of the Horndeski type via the Einstein tensor. 
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6 Dark Matter-Dark Energy non-minimal couplings? 

In this section we will discuss another intriguing possibility of NMC fluids. In fact, one can 
speculate whether there could be more than one fluid which is NMC. In particular, one can 
imagine the situation in which there are two NMC fluids one playing the role of DM and 
another one related to DE, such that couplings of the form 

Fi{ni,si)u^R^yv!{ + F2{n2, S2)u'^R^uU2 + Ei2(rei, si; n2, 52)^1 > (6-1) 

can be available. The last term in particular, represents an interaction, mediated by gravity, 
between DE and DM fluids. Recently, it has been shown that Planck data allow for this 
possibility [40] or even more, they seems to favour such interaction [41, 42] so that coupling 
as those presented above may represent an interesting extension to the model presented in this 
paper able to include such experimental hints. Interactions of this kind are also interesting in 
connection to other models like the one investigated in [43] , where a Lorentz breaking vector 
field is coupled to a fluid DM fluid, or the one investigated in [44] where a vector version of 
the Horndeski action is constructed. 

If we consider the simplest model in which only the NMC between the two fluids is 
present we get the Einstein equations 


dF, 

, I n. _n _ /y’, 1 _U I nr 




+ {R) + hW T hg) 


dni 


n2 


dm 




= 0 , ( 6 . 2 ) 


where rgj is the SET for the i-th component as given by (2.3), while 

= (Frfttgug)/2. The fluid are better derived from the direct variation of the action 
with respect to the fluid variables of both components. In particular, unless some mixing 
between the constraint is introduced, both entropies and number densities will be separately 
conserved. On the other hand, the Euler equation gets modified as 




5Fi 


12 


dm 


Fi2 

m 


+ h7y 


— —ai2h' 




(*) 


(*) "P 

V, 




5Fi 


12 


-T12 
m ' 


iliaUj 


dm 

- V, 


_ 

m 
Fi2 


m. 


RpivU^j 


(6.3) 


where i,j = 1, 2, i / j and F* = Fi{m) and F 12 = Fi 2 (ni, 77 - 2 ). 

To summarize, in the case of two interacting fluids, we end up with a system in which the 
two fluids have conserved charges along the flow, m and Si and a modified Euler equations. 
In particular, the inhomogeneities in each of the fluid act as force terms for the other. 


7 Discussion and conclnsions 

In the era of precision cosmology and with many upcoming cosmological and astrophysical 
surveys, alternatives to the ACDMwill be dramatically put under scrutiny and most likely 
we will soon be able to distinguish among the different proposals. In this spirit and motivated 
by the quest for a better understanding of the dark component of the Universe and of the 
fluid dynamics in curved space-times, we have investigated the full dynamics of a NMC fluid 
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where both a conformal and a disformal coupling to curvature are present. We have seen 
how the requirement of no fluid derivative in the action reduces the possible NMC only to 
two terms, FcR and ■ These two terms are however general enough to produce 

significant modifications with respect to the simple minimally coupled case. 

In particular, the Einstein equations are much richer and shows the appearance of 
derivative terms of fluid variables which play an important role in the Newtonian limit. 
In fact, the Poisson equation may be sourced by gradients of the fluid density distribution 
meaning that in these scenarios it will depend not only on how much matter is there at a 
point but also on how it is distributed. We have also pointed out how the conformal and the 
disformal couplings have different features: while the former is related to anisotropic stresses 
the latter is not, thus it does not excite new gravitational degrees of freedom with respect to 
GR. These features directly relate to observable quantities and hence are potential tools for 
constraining this class of models. 

For what concerns the fluid equations, we show how in the Euler equation appears an 
extra force term related to the curvature, which is a clear and distinguishable feature of 
the NMC, while the continuity equation is modihed in such a way that it will still allow 
the definition of conserved quantities along the flow. Hence, structure formation will be 
affected by the NMC in a way that may potentially stabilize their growth under a certain 
scale (as a consequence of the extra force). On the other hand, as one can see from the 
Friedmann equations, changes are expected to occur also at the background level. The 
model can be easily made compatible with early universe physics, which is well described by 
the ACDM model, while still providing interesting new features at present. 

We have also discussed the important issue of extra degrees of freedom that may be 
propagated in this kind of models as a consequence of the NMC. In section 3.3 we have 
provided an argument for why this should not be the case. Of course, only a complex Hamil¬ 
tonian analysis can conhrm this conclusion but we reserve this analysis for a future work. 
It is also important to stress that the model presented in this paper is a phenomenological 
one and not a fundamental theory of gravity and that it has been though as a cosmological 
extension to the ACDM model. For this reason, in section 4.2, we have investigated the dy¬ 
namics of linear perturbations around a cosmological background and showed how it is free 
of extra degrees of freedom. We hence conclude that in the regime for which the model has 
been thought to apply to there are no unwanted degrees of freedom. 

Finally we have considered the case of two NMC fluids which can interact via curvature 
couplings. This is an intriguing possibility as it may allow DM-DE interactions an interesting 
extension that seems to be favoured by current cosmological data. 

The issue of non-minimal couplings in gravitation theories is a long standing one. When 
dealing with cosmological fluids we have argued that non-minimal coupling for exotic com¬ 
ponents like DM and DE should not and cannot be neglected (even just to be able to exclude 
them on the base of futnre observations). We have presented here an exploration of the 
implications of this extension of the standard cosmological framework and found several in¬ 
teresting novelties that we think should stimulate deeper investigations of these models. We 
think that NMC fluid models are now mature for phenomenological applications and we hope 
to further advance in this direction in future. 
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A Equation of motion from the fluid variables 

In this appendix we report the equations derived from the NMC actions (3.6) and (3.19) when 
the variation is taken with respect to the fluid variables. This has been used as a consistency 
check that the equations obtained from the covariant conservation of the metric equation 
and those derived from a direct variation of the action with respect to the fluid variables 
are equivalent. Moreover, the equations presented in this appendix help in clarifying some 
thermodynamic aspects of the NMC fluid, in particular making explicit the way they modify 
the standard relations. 


A.l The case of conformally coupled fluid 

The variation of the action (3.6) with respect to the fluid variables gives the following set of 
equations 


0 

0 

0 

0 

0 

0 


5S 

6J^ 

d(p 

Jo' 


fdF dFA 

\on on J 

-V^(sJ^), 

, fdF BFc 


6S 

5S 




+ V^(/9 + , 


+ , 


(A.l) 

(A.2) 

(A.3) 

(A.4) 

(A.5) 

(A.6) 


Only two equations are modified with respect to those obtained from a minimally coupled 
fluid (2.9)-(2.14) and as it has been discussed the modiflcations occurs as rescalings or shifts 
of existing fluid variables. Eor example, the flrst equation, contracted with the four velocity 
u^, and the fourth one, reads 


- (F + OcRFc) = ^p + s9 , 
on 

--^{F + acRF,) = 9. 


(A.7) 

(A.8) 
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Given the similarity if the right hand side of these equations with those for the minimally 
coupled case it is straightforward to make the following thermodynamic identifications 

Pc = —{F + acRFc)i (^-9) 

V9 = /, (A.IO) 

e = T. (A.ii) 

where T is the temperature and / is the chemical free energy, as defined in section 2. 

From the set of equations (A.1)-(A.6) it is easy to show that one gets the same equations 
as from the conservation of the Einstein equations. In fact, from equations (A.2) and (A.3) 
we recover equation (3.12), while from equation (A.l), using the procedure described at the 
end of section 2, we recover exactly equation (3.15). Hence, we conclude that the addition 
of the conformal NMC is compatible with the assumption to have a perfect fluid. 


A.2 The case of disformally coupled fluid 

We proceed now with the computation of the disformally coupled equations from the action 
(3.19). The variation with respect to the fluid variables gives the following system of equations 
is 

(5 S dF 2 / dF \ 

0 = — = + sV„e + - -Od ( “ FdR^yu'^h^^A , 


(A.12) 

0 = ^ = -V^J^ (A.13) 

0 = ^ = -V4sJ^), (A.14) 

0 = ^ = -V4/3^J^), (A.16) 

0 = ^ = J^V^a^, (A.17) 


where as for the previous case of conformally couple fluid only two equations are modified. 
These modifications occurs in a very similar way as what happens in the case of the confor¬ 
mally coupled fluids, the only difference being the last term in the first equations. However, 
when contracted with the four velocity of fluid such term disappears. Hence, we can define 
as in the previous case the following thermodynamic quantities 

Pd = ~{F + 2ad{R)Fd) i (A.18) 

<p = f , (A.19) 

0 = T. (A.20) 

where again T is the temperature and / is the chemical free energy, as defined in section 2. 

Also, using the same procedure described above one can show that the equations for the 
fluid (A.12)-(A.17) implies those obtained form the conservation of the Einstein equations. 
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B Linear perturbation coefficients 


In this appendix we report the detailed expressions of the linear perturnation coefficients 
used in section 4.2. 


B.l Conformally coupled fluid 

The coefficients of the space space component of the Einstein equations are 
^si — {3H^ + 2H) + ac(^2,H[2pc + 3pc) + H‘^(3pc — 9c^(pc + Pc))^ > 

Es 2 = H + acH{2pc + 3pc), 

Es3 — (3i?^ + 2H + K^p) + OLc{2H(2pc + 3pc) + H‘^{3pc — 9c'^{pc + Pc))^ j 

Esi = H - acH{-pc + 3cl{pc + Pc)) , 

Es 5 = 1 + ac{2pc + 3pc) , 

Esq = [3H - (1 - 3cl)H) {pc + Pc), 

Es7 = H^{2 - lOc^ + 9c^) + 3H{1 + 2cl). 


The coefficients in the Euler equation are 

Esci = P E 3ac{pc E Pc){2H'^ + H), 

Esc 2 = P E 3ac{pc E Pc) I 2H^{1 — 3c^) E H(5 — 3c^) + — 


Esc 3 = P- 3aci2H^ + H){pc + Pc), 

Esc4 = {2H^ E H ). 


B.2 Disformally coupled fluid 

The coefficients for the Euler equation read 

Edi = pEp-ad{-H{pd-3pd) E 3H‘^ {pd E Pd)) , 

/ ii 

Ed2 = (1 - 3cl){pEp) - ad(^- — {pd - 3pd) 

— 3(—1 + 3c^)H^{pd E Pd) E — 9c^)pd E 3(7 — 3c'^)pd)^ , 
Ed3 = pEp-ad{H{pd - 3pd) - 3H^{pd E 5pd)) 

Coeffcienents of the space space component of the Einstein equations 

^sdl — 2 “1“ 3pd) ) 

Esd2 = + 4ff + Od (3cdPd - (2 - 3c\)pd) - 2H[pd E 3pd)'^ , 

^sd3 — ^ds2 1 

EsdA = 6H - 3adH[2pd - 3{pd EPd)c^d) > 

Esdb = 2H - adH{pd E 3pd ), 

Esdd 1 '^{Pd E Pd) ) 

Esdi = H[2pd - 3(1 - cl){pd Epd)) , 

Esds = ^dE H {l E 3(?d)) 


(B.l) 

(B.2) 

(B.3) 

(B.4) 

(B.5) 

(B.6) 

(B.7) 

(B.8) 

(B.9) 

(B.IO) 

(B.ll) 


(B.12) 

(B.13) 

(B.14) 

(B.15) 

(B.16) 

(B.17) 

(B.18) 

(B.19) 

(B.20) 

(B.21) 

(B.22) 
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